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We present a supersymmetric D-brane model that has CP spontaneously broken by discrete 
torsion. The low energy physics is largely independent of the compactification scheme and the 

. kahler metric has 'texture zeros' dictated by the the choice of discrete torsion. This motivates 

a simple ansatz for the kahler metric which results in a CKM matrix given in terms of two free 
parameters, hence we predict a single mixing angle and the CKM phase. The CKM phase is 

. predicted to be close to |. 

m : 

O ■ 

1 Introduction and motivation 

<n : 

CP violation is a curious aspect of beyond the Standard Model physics. Baryogenesis (or rather 
its absence in the Standard Model) almost certainly indicates additional sources of CP violation 
beyond that which has been established by experiment to exist in the CKM matrix [I], p|, [|. On the 
other hand the absence of electric dipole moments (EDMs) seems to indicate that the additional 
CP violation must have a very constrained form, possibly dictated by a direct connection with 
flavour structure. 

^ String theory may be giving us an important clue as to the real nature of CP violation; in string 

theory, what we call 4 dimensional CP is actually a gauge transformation plus Lorentz rotation 
of the 10 dimensional theory ||, |5|. Thus string theory is one of a finite class of theories in which 
CP is a discrete gauge symmetry which can only be spontaneously broken ||. Encouraged by 
this observation, a number of authors have attempted spontaneously to break CP in the effective 
supergravity approximation to various string models || [7[ §, ||, [l(| [TT], [12], |13|]. However, in these 



models it has proven very difficult to find a satisfactory suppression of EDMs and flavour changing 
processes. The problem seems to be that whatever fields break CP and generate flavour structure 
also contribute to supersymmetry breaking (specifically, their auxilliary fields aquire non-zero F- 
terms). Thus the supersymmetry breaking 'knows about' the flavour and CP structure, giving 
rise to the usual supersymmetric flavour and CP problems. 

We think that these problems can be avoided if one instead establishes the flavour and CP 
structure at the string theory level rather than the supergravity level. This is because string theory 
allows us to maintain supersymmetry which should in turn allow us to separate CP and flavour 
from supersymmetry breaking. In this paper we take a first step in this direction by constructing 
a supersymmetric MSSM-like model which has a full flavour structure and broken CP. In a later 
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paper we shall examine the consequences for supersymmetry breaking, although we shall make 
some preliminary observations here. 



Our framework will be the 'bottom-up' approach to string model building developed in ref [JT4 
This approach allows one to construct a local configuration of D-branes that reproduces most of 
the phenomenological features of the MSSM without having to worry about the global properties 
of the compactification. The source of spontaneous CP violation will be discrete torsion, a choice 
of orbifold group action on the B-field background, analogous to Wilson lines for gauge fields ||T5|j . 
One nice feature of discrete torsion is that the Yukawa couplings contain a non-trivial complex 
phase, with the breaking of 4 dimensional CP appearing as a phase in the CKM matrix. With 
a rather simple ansatz dictated by our choice of discrete torsion (similar to 'texture zero' models 
in the MSSM) we find a CKM matrix defined by two free parameters. As a result, we predict 
a single mixing angle and the CKM phase, both of which lie within current experimental limits. 
Furthermore, quark mass ratios are determined and are seen to be reasonably close to experimental 
values. 

It turns out that the up-quark Yukawa couplings are hermitian. The general idea of hermiticity 
has often been proposed as a way of suppressing EDMs in softly broken supersymmetric models |16 



17[ fL9| , |20H but it has been quite difficult to achieve in field theory. Here it arises as a natural 
consequence of the form of the Yukawa couplings with discrete torsion, indicating a possible 
resolution of the supersymmetric CP problem. 

The paper is organized as follows. In the following section we introduce the model which is a 



generalization of the discrete torsion models presented in ref [|14|]. Particular attention is paid to 
the effect of discrete torsion on the superpotential and it's role in breaking 4 dimensional CP. We 
then proceed to a more phenomenological discussion, including Yukawa couplings and the CKM 
matrix. 

2 A D-Brane model with Discrete Torsion 

First, a brief descriptive overview is given, and then a more detailed account that concentrates on 
the form of the superpotential. 

2.1 A brief overview 

The model we are interested in contains the following ingredients, 

• A non-compact internal space C 3 , containing a Z 3 x Zm x Zm orbifold singularity. 

• Coincident D3-Branes whose world- volume theory has the gauge group SU(3)c x SU(2)l x 

U(l)y. 

• Mutually perpendicular D7-BranesQ to cancel tadpoles, each with gauge group U(l) x U(2). 

The D3-branes and D7-Branes are located on the Z 3 x Zm x Zm singularity, this breaks 
supersymmetry from N=4 down to N=l, allowing for chirality. The orbifold group leads to model 

1 We call a D7-brane whose worldvolume is Zi — a D7i-brane. 
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Figure 1: The D3-branes are shown separated for clarity, although they are in fact coincident, 
and the D7-branes fill the figure. The position of end points of the strings determine how they 
transform under the SM gauge group, with fields such as the right-handed leptons and down 
quarks having one end on a D3 brane and another on a D7-brane. 

with three generations of Higgs' in addition to the standard quark-lepton generations. Furthermore 
the Ijm xZ tf factor allows for a choice of discrete torsion (this choice is labelled by an integer s). 
The discrete torsion then gives rise to a complex phase in the superpotential, these facts will be 
discussed in more detail in the rest of the paper. 

Closed strings propagate only in the bulk and are gauge singlets, whereas the open strings are 
localized to the D-branes. The D3-branes are arranged in sets of 3s, 2s and s corresponding to the 
gauge groups U(3), U(2) and U(1)R. Open strings with ends on D-branes within the same set give 
rise to the gauge bosons, and those with endpoints in different sets of D-branes give rise to matter 
fields. This makes it very simple to identify the MSSM fields as illustrated in figure [L| The gauge 
group on the D3-branes is £7(3) x £7(2) x £7(1), remarkably however, only one linear combination of 
the three U(l) factors is non-anomalous and this combination gives rise to the correct hypercharge 
assignments. This results in a world-volume theory with standard model gauge group. 

2.2 Orbifold action on the string spectrum 

A system of overlapping D3-branes and D7-branes has the following massless open string states 
in the spectrum, 

1. In the 'untwisted' 33 sector, 

• world-volume gauge bosons A 33 ^^; 1//2 10 >, 
2 The action of Zm x Zm on the string spectrum breaks the gauge group U{sn) down to U(n). 
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• complex scalars[] A33 , 0!_ 1 / 2 |0 >, i=l ,2,3. 

• Weyl fermions A 33 |si, s 2 , s 3 , s 4 >, Sj = ±1/2 and J2t=i s i °dd. With s 4 determining the 
space-time chirality. 

2. In the 'twisted' 37j sector, 

• complex scalars A 37i |s,j, >, j, ^ {^,4} and = s&, 

• Weyl fermions A 37i |sj,S4 >, Sj = S4. 

Here Ajj are the Chan-Paton matrices. 

The action of the orbifold on these states can be split up into it's action on the Chan-Paton 
matrices and the action on the string worldsheet fields . For aZjv orbifold with generator 9 and 
twist vector Ve = (ai, a 2 , a 3 )/iV, N=l supersymmetry requires a\ + a 2 + a 3 = 0. The generator 
9 acts on the 33 sector by leaving the gauge boson invariant, and transforming the scalars and 
fermions respectively as, 

^_i/ 2 |0> = e 2 ™^ 1/2 |0>, (1) 
9 k \s> = e 2mg - vl >\s>. (2) 

Here s — (s 1; s 2 ,s 3 ), and the twisted 37j states transform as the fermions above if we set any Sj 
not defined to be zero. For our particular case we take the twist vectors for Z 3 and the two Zm 
factors to be, 

v e = (l,l,-2)/3, 
v m = (1,-1, 0)/M, 

v U2 = (0,1,-1)/M, (3) 

respectively. The action of 9 on the Chan-Paton matrices of the twisted and untwisted states 
respectively is, 

9\ 3 3 = 7e,3A 33 7e j3 , (4) 

0A 37 . = 70,3A 37l 7j, v (5) 

where 

7e i3 = diag(I 3s ,aI 2s ,a 2 I s ), (6) 

70 i7i = diag(al s ,a 2 l 2s )- (7) 

Here we have defined a = e~ . The 70 matrices form a unitary representation of Z 3 determined 
by our requirement of a SM gauge group and the tadpole cancelation condition, 

3 

^Tr(70 )7i )+3Tr(70, 3 ) = O. (8) 

i=l 



3 These are the collective coordinates describing the embedding of the D3-branes in the internal space. 
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The action of the generators Ui is the same as above, however this time the j Ui matrices form 
a projective representation of Zm x Za{ |^T|, ^], a possibility that arises since the action of 



Zm x Zm allows for a choice of discrete torsion. This is defined in terms of a cocycle (3(g, h) e 



x Z M ,U(1)) S Z 



p3j |24| , which manifests itself in the projective representation where 
Iglh = fl(g, h) r Yh'J g - It follows that there are M non-equivalent choices of discrete torsion, described 
by distinct cocycles (3(u>i, where n=l,..,M. Such possibilities are conventionally 

distinguished by an integer s = M/gcd(n, M). For a string endpoint residing on a generic set of 
sn W D-branes, the action on the Chan-Paton factors is determined by the matrices, 



7wi = ®l,mU M %i ® 
7^2 = ®l,mUJ M % 2 



7}.: ' ' 



(9) 
(10) 



where u M = e**' , m n}£ = n®, I, m = 0, . . . , (^) - 1, and 
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(ii) 



(12) 



Since Tr(7^.) = tadpoles are unaffected by discrete torsion, it follows that there are no restric- 



tions on the n[^. For simplicity, we choose, 



n 



W 

lm 



n® if l=m=0, 
otherwise. 



(13) 



Combining the above transformations results in the projection equations which determine the 
spectrum and gauge group of the model to be that discussed above. For example, in the 33 sector 
we have, 



A33 — 76»,3^337e,3 
A 3 3^ii /2 |0 > 

A 33 = e 27ri / 3 7 (?i3 A337 ( !,3 

A33^!i/ 2 |0 > 
^33 — 

A 33 V' 3 1 -I0> 



A33 — 70)1,3^337^1,3 



A33 = 7o) 2 ,3A 3 37 



^2,3 



A33 = e ^7o)i,3A337ii,3 A33 = 7 W2 ,3A337 [ 



t 

u>2,3 



(14) 



e 2W3 7e,3A 3 37e,3 A 33 = e Te^sAssT^ A33 = ^ Twa^AasT, 
In • 

1/2 



A33 = e 47n/3 7e,3A 3 37j, 3 A 33 = Tui.aAaaTli.a A 33 
Similar projection equations hold for superpartners and the 37« sector. 



^2,3 



A33 = e u 7o) 2> 3A 33 7i 2 ,3- 
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2.3 The shift formalism and the spectrum 



A particularly useful formalism for finding the multiplets is the shift formalism of ref.[|2|, [25|, |26], £7| 
which makes use of the fact that the basis we have chosen for the gamma's is the Cartan basis (i.e. 
it is the basis where the diagonal elements correspond to the generators of the Cartan subalgebra). 
We write the gamma embeddings more explicitly^ 

7e , 3 = diag[(l) sn \(a) sn \(a 2 ) sn >] (15) 
7*1,3 = diag[(l, e~\ ...e'^T , (1, e~\ ...£-^) n \ (1, e" 1 , ...e~^) n '] (16) 
7^2 ,3 = blockdiag [ 7w2 In > 7^2 -^ni ? 7^2 -^«2 

An obvious similar structure holds for the 77's which we will not show here. Now the diagonal 7 
matrices we write in terms of shift vectors that act on the root vectors of the states p; 

7*,3 = ^ (18) 
7.1,3 = e™ v ^- H > (19) 

whereas (for this choice of torsion) the effect of 7^ is a pure permutation of the root elements. 
For the gauge states the projections are (mod 1) 

p.Ve = (20) 
p.V Ul = (21) 
IL 2 (p) - p, (22) 

where the final equation means that we must build a linear combination of the remaining states 
which is mapped into itself under permutations. The following linear combinations of roots satisifes 
eq.© 

p= [ ((±0..0)(=F-0)(0-0) w °~ 2 )(0-0) Wl (0-0) W2 ] (23) 

+ [((0±0..0)(0^0..0)(0..0)"°- 2 )(0..0) ni (0..0) n2 ] + ... (24) 



p= [ (0..0) w °((±0..0)(T..0)(0..0) ni ~ 2 )(0..0) n2 ] (25) 

+ [ (Q..0)"°((0 ± O..Q)(0 =f 0..0)(0..0) ni " 2 )(0-0)" 2 ] + ■■■ (26) 

p = [ (0..0) no (0..0) ni ((±0..0)(T--0)(0-0) n2 ~ 2 )] (27) 

+ [ (0..0)"°(0..0) ni ((0 ± O..Q)(0 =f Q..0)(0..0) n2 ~ 2 )] + ... (28) 

where the underlining means all permutations of brackets. Once we add the Cartan generators 
we find gauge fields in the adjoint representation of U(no) x U(n\) x [/(712). For the chiral states 
the projections are (mod 1) 

P-Ve = v 0! i 

n W2 (p) - e 2 ™»^p. 

4 Previously we chose no = 3, ri\ = 2 and n-i — 1 to obtain the SM gauge group but here we leave them general. 
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The multiplets are all in the bifundamental in this case and are 



(r^,ni,l) + (l,nT,n 2 ) + (n ,l,r^) = (3, 2, 1) + (1, 2, 1) + (3, 1, 1); (29) 

i.e. with n = 3, n\ = 2 and n 2 = 1, we find a left-handed quark doublet, a right-handed quark 
singlet and a Higgs doublet. As an example consider bifundamental (no, n\, 1) states coming from 
the 3 rd complex plane, i = 3. To satisfy the projection equations the roots must be the following 
linear combination 

p = [ (-o..o)(o..o) ra °- 1 (+o..o)(o..o) ni - 1 (0-0) n2 ] (30) 

+ e [(0 - 0..0)(0..0) w °- 1 (0 + O-.O)^.^)" 1 - 1 ^.^)" 2 ] (31) 
+ e 2 [ (00 - 0..0)(0..0) n °- 1 (00 + 0..0)(0..0) ni - 1 (0..0)" 2 ] +... (32) 

Because the projection equations for the u W2) j are the same for all i, we find the same multiplets 
from each complex plane, and so for these models we have three generations. 



2.4 The superpotential 

The form of the superpotential can be deduced as in |28| and includes the following terms, 

2 

W = e abc Tr(<$>%$ b 3 <$> c 3 ) + ^Tr($ 3 $ 37l $ 37 *) + . . . , (33) 

i=0 

where $3 = and $ 37i = A 37i y? 371 , with <p 3 a 33 chiral superfield and y? 37i a 37j chiral 

superfield. Using the expressions derived for the Chan-Paton factors in appendix [A], we have, 



= e^Tr(X& ) X<gxg)Tr(Z a Z b Z°) + ..., 



where Z % = Y^(p 3 . Redefining the VW's by a constant such that a factor of Tr(X 3 ^ X 3 ^ X 3 f) is 
absorbed into each term and noting that 

we arrive at, 

W = e abc Tr(Z a Z b Z c ) + (1 - e - 27Ti/M )s abc Tr(Z a Z b Z c ) + ... (36) 

Here s 132 = s 321 = s 213 = 1 and all other components are zero. Identifying our supersymmetric 
SM fields, equations ( |33"|) and ( |36| ) give us the yukawa type couplings, 

W = h abc TF R Q b L H c u + f abc T R Q b L H c d + . . . , (37) 

where h abc = e abc + (1 — e~ 2m l M )s abc and f abc = S ab 5 bc . Hence we have an up-quark Yukawa term 
containing a complex phase which gives rise to CP violation, and a prediction for the phase of the 
CKM matrix is presented in detail in section |||. Before deriving the CKM matrix however, let us 
discuss the nature of the CP violation. 



2.5 More on CP violation 



In the heterotic string, it is well known that a four dimensional CP conjugation corresponds to a 
rotation in the full 10 dimensional theory, and thus 4 dimensional CP turns out to be a discrete 
gauge symmetry |4], || (i.e. part of the 10D Poincare group). The situation is heuristically as 
shown in fig[2|. The three non-compact space dimensions are labelled Xi, with yi and Wi labelling 
6 internal compact dimensions. Parity is defined as a reflection in one direction, X\ say. By a 
rotation through n in the X2X3 plane, four dimensional parity is equivalent to a reflection in all the 
space directions, the more conventional definition (valid when the number of space dimensions is 
odd). However if we simultaneously reflect in an odd number of internal compact dimensions the 
total is a rotation in three orthogonal planes. Complexifying the internal space as 

Zi = Wi + iyi, (38) 

we can, without loss of generality, choose the additional reflection corresponding to a CP trans- 
formation to act on the i/i coordinates, with 

CP : x» -> x M , Zi -> z*. (39) 

By supersymmetry we must also reflect the fermionic superpartners, so that the transformation 
reverses the R-charges of the spectrum. In turn, the requirement of modular invariance is only 
satisfied if we reverse the gauge charges as well. (In the standard embedding this is a simple 
copying of the internal reflections on the space-time side into the SU (3) subgroup on the gauge 
side.) The nett result of the combined rotation plus gauge transformation is a 4 dimensional CP 
transformation. 

An even simpler picture holds for branes at orbifold fixed points, as we shall now see, and it 
is the geometric configuration of torsion plus singularity that breaks CP. As an example, consider 
the chiral multiplets written explicitly in eq. ([30|) . The corresponding antiparticles are given by 
the conjugate excitation, and the projections are reversed accordingly; 

-v ,i (40) 

-*W ( 41 ) 
e~ 2mv ^% (42) 

giving 

p = [ (+0..0)(0..0)" - 1 (-0..0)(0..0) Wl - 1 (0..0) n2 ] (43) 

+ e- ^O + 0..0)(0..0) w °- 1 (0 - Q..0)(0..0) rai - 1 (0-0) n2 ] (44) 
+ £~ 2 [ (00 + 0..0)(0..0) n °- 1 (00 - 0..0)(0..0) W1 - 1 (0..0) W2 ] +... (45) 

as expected. 

Now consider the rotation x M — > x^, Zi — > z*. Again by supersymmetry we transform ipi — > ipi 
and need to reverse the chiralities in the Ramond sector. The particle projections now trivially 
become antiparticle projections. To get the right antiparticles however, we also need to ensure that 
we have the same projective representation which is the case if the discrete torsion is unchanged 



P-Ve = 

P-Van = 
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Figure 2: Rotations on 3 non-compact plus 3 compact directions corresponding to CP. 



(e.g. so that we have same linear combination as in eq.fl43|). If this were not the case we would 
simply be unable to identify these particular rotations with CP conjugation.) For open string 
sectors the effect of the torsion in a sector with Zm x Zm twists given by (g = a, b ; h = a', b') 
can be written as 

f3(g,h) = e ab '- ba '. (46) 

The effect of the z,- t — > z* reflection on the torsion in any sector is simply to reverse the twists, 
(a, b ; a 1 , b' ) — > — (a, b ; a', b' ) so (3(g, h) is indeed unchanged and the rotation gives us the antipar- 
ticle. i.e. 4 dimensional CP conjugation does indeed correspond to a rotation in the 10D theory, 
with or without discrete torsion. 

It is precisely because the discrete torsion is invariant under CP conjugation that the super- 
potential interactions are not. Indeed the superpotential is of the form 

W = Tr($i$ 2 $3 + £~ 1/n $ 2 $i$3), (47) 

where here the $'s stand for generic chiral fields, and where a summation over the SU(3) Levi- 
Cevita symbol is implied. Under the CP transformation we conjugate the chiral superfields but 
the discrete torsion remains the same. Thus 

J ' dd 2 W + J dfW -> J S 2 W C p + J d6 2 W C P (48) 

where 

W CP = Tr($!$ 2 $ 3 + £ 1/n $ 2 $i$ 3 ). (49) 

So it is the discrete torsion term in the superpotential that is directly responsible for breaking 4 
dimensional CP and, because CP conjugation is associated with a 10D Lorentz rotation, it is a 
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spontaneously broken symmetry. We will now divert our attention to the more phenomenological 
issues. 

3 Yukawas and the CKM matrix 

In this section we first discuss some theoretical considerations for obtaining the canonical yukawa 
couplings. In particular, we derive bottom-up properties of the kahler metric that have not been 
previously noted. We then proceed to the CKM matrix and predictions for the CKM angles and 
quark mass ratios. Finally we discuss radiative corrections between the string and weak scales. 

3.1 Normalisation of Yukawas 

3.1.1 Bottom-up properties of the kahler metric 

To be consistent with N=l supersymmetry, we assume the low-energy effective field theory de- 
cribing our visible sector is an N=l SUGRA theory. It is well known that in such a theory scalar 
fields have a non-canonical kinetic energy term given by a kahler metric which is a function of the 
compactification moduli. This normalistion of scalar fields, results in an effective normalisation of 
coupling constants. In particular, if we expand the kahler potential to first order in the charged 
matter fields^, we have 

K = K(M a , M a *) + K(M a , M a *) a/3 Tr($ a *$ /3 ) + (50) 
where M a , M a * are the closed string moduli fields and we obtain, 

kbc = e* /2 /i, mn (K~ 1/2 ) , a(K" 1/2 )' n 6 (K~ 1/2 ) n c , (51) 

where h abc are the canonically normalised yukawa couplings. It is important to note that K Q/ g 
factorises into different components, one for each sector of the spectrum, e.g. in our case the 33 
and 37 sector fields have the seperate kahler matter metrics, (K 33 ) aj g and (K^)^. 

One of the strengths of the bottom-up approach is that, up until now, our discussion has been 
independent of any compactification scheme. In particular, the gauge group, SM spectrum and 3 
generations of matter fields are all features independent of the global structure of the model. This 
affords us two approaches to the determination of the yukawas, 

• Select a particular compactification and calculate the exact kahler potential. 

• Choose an ansatz for K 33 and K 37 (note that the CKM matrix is independent of K). 

The first approach is technically difficult and only acheivable in the simplest of cases, e.g. toroidal 
orbifolds. In addition, it goes against the bottom-up approach, which is to constrain the model 
as much as possible from the local rather than global properties. The second approach, at first 
sight, seems too arbitrary. However, the choice of possible ansatze is restricted by two observations. 
Firstly, the discrete torsion can constrain the form of K 33 . Such a situation arises as a consequence 

5 i.e. fields charged under our visible sector gauge group. 
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of Tr($^$f 3 ) containing a factor Tr{{X$)iX$) arising from the Chan-Paton factors. This 
matrix contains 'texture zeros', dependent on the discrete torsion(n), which are inherited by 
(K33)q,/3. Note that the kahler metric is hermitian, and we will coventionally refer only to the zeros 
in the upper triangular section of K33. There are three possible cases, 

• minimal discrete torsion, i.e. n = 1, with K33 diagonal, 

• n = 2, with two off diagonal zeros, 

• and n > 3, with no off diagonal zeros. 

These results are derived in the appendices. 

Finally, we have the restriction that K37 is diagonal. (K 37 )ij can be calculated from tree-level 
modulus-matter scattering [p9|1 . The amplitude is given by a disc diagram with two open string 
vertex operators on the boundary, a (37«) and (7^3) vertex, and two closed string vertex operators 
attached to the interior, corresponding to the moduli fields. However vertices of (37j) and (7j3) 
strings must come in pairs in order to get a consistent D-brane boundary on the disc. It follows 
that in this case i=j and hence K37 is diagonal. 

This allows us to explore possibilities for interesting flavour structure and CP violation in our 
yukawas without worrying about the exact compactification scheme we are using. That is, we 
choose an ansatz for the kahler metric based on our choice of discrete torsion, not our choice of 
compactification. 

3.1.2 Explicit expressions for Yukawas 

Defining the hermitian matrices, 

t33 = K3 V (52) 
t 37 = K 37 , 

and substituting into (|5ll) we obtain, in the general case, canonical up-quark yukawas, 

Uabc _ KI2-L (A/2\la(4l/ 2 \mb(A/ 2 \nc 

IL — C lHmn\ L 33 ) ^33 ) ^33 ) /rn\ 
_ p K/2( r , (1 _ „-2m/M\ a \(A/2\la(A/2\ m b(A/2\nc V JO / 

— e yt-lmn 1 \ L e J*/mnA''33 / \ L 33 I \ L 33 I 

Furthermore giving the three generation of Higgs fields vevs, (h u ) c , we obtain 

(Y u y b = e k l\e lmn + (1 - e- 2m / A ^ Sim 0(^ 2 )' a (^ 2 ) mb (43 2 ) nc (^)c. (54) 



This expression also holds if we associate the actual physical Higgs fields with a linear combination 
of the Hl Jd fields. That is 

K /d = K%Hi /d , (55) 

where H^u are the physical Higgs fields and U u /d is a unitary matrix. Then assuming that only 

Hh d is light and the other two physical Higgs' are heavy, we simply replace (h u ) c with (U^) c i- 
Similarly, for the down-quark yukawas, we obtain, 

(^"^e^X:^ 2 )^^ 2 )^^ 2 )^^), (56) 
i=i 
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Using these formulae we can calculate the CKM matrix and determine the possiblities for CP 
violation in our model. 



3.2 The CKM matrix 

3.2.1 A simple ansatz from a choice of discrete torsion 

Using the bottom-up properties of our kahler metric, that is, 

• the zeros in K33 are determined by the discrete torsion, 

• K37 is diagonal, 

we try the following simple ansatz, take M=5, n=2 and 




& = 1 0, (57) 



q 7 = diag[l,£,l], (58) 

ih )= fr 1 ' }' (59) 
(h d ) = (l,e,e). v ' 

Notice that our choice of discrete torsion introduces two off diagonal zeros in K 33 which are 

1 

preserved in t| 3 . 

Using this ansatz we can calculate the Yukawas and CKM matrix. The Yukawas have the 
following form, 

/ (l-e- 2i7T/5 )r]e 2 (1 - e -2 ^/ 5 )^ 3 (1 - r]e~ 2i7T/5 )e 
Y u = a 2 e (1 - e- 2i ^ 5 )r]e 3 (77 - e~ 2i7T / 5 )e 2 | , (60) 

\ (rj- e~ 2 ™/ 5 )e (1 - r]e~ 2i7r / 5 )e 2 (1 - e ~ 2l7T / 5 ) 

and 




Y d = as f- . (61) 



With the freedom to rephase the right handed quark fields (because the model is supersymmetric 
we can make rephasings just as in the Standard Model), Y u can be put into hermitian form. Notice 
that a appears only as an overall factor in Y u and in a single entry in Y d , as a consequence the 
CKM matrix is only dependent on the two parameters e and r], with a determining the mass of 



the down quark. Up to overall factors, which include the prefactors of (|60| ) and (|6l|), the masses 
of the quarks are given approximately for small r\ and e by, 



e 



2 



m u = l.9 V e\ mc = ^le 2 -^ V e 2 + ^l V 2 e 2 , m t = ^ + jje 2 - ^ V e 2 

The CKM matrix can be written approximately in terms of these mass eigenvalues, but it is 
not very illuminating to show the functional dependence here. Instead, in the next subsection, we 
depict graphically it's functional dependence on e and rj. 
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3.2.2 Predictions: CKM angles and Quark mass ratios 



In the standard parametrisation, the CKM matrix is given in terms of three mixing angles, sin # 12 , 
sin #23, sin0 13 , and a complex phase 5. Figures 3 to 6 illustrate the dependence of these CKM 
angles on e and r\. The blue shaded areas highlight the regions where the CKM angles agree with 
experiment |3(|. Since we are discussing string scale predictions we must also include RGE effects 
in running from the string to the weak scale. This effect is discussed in the next subsection, where 
we estimate that sin #23 and sin #13 are suppressed by approximately 30%. This is depicted in the 
figures by a red shaded region which higlights values of e and r\ that will reproduce experimental 
values after taking into account RGE effects. This information is summarised in figure ^ where the 
regions in parameter space which produce a realistic CKM matrix are highlighted. In particular, 
it can be seen that S is approximately tt/3. 

The quark mass ratios are also determined, ^ can be fixed by choosing a value for a, while 
the rest of the quark mass ratios are dependent only on e and 77. Average values for the blue 
and red regions are presented in table [I] and compared to experimental values. As can be seen, 
we have good agreement with experiment except for ^ which is too large. Note that we have 
not included RGE effects here, which would tend to increase the value of — . As we have chosen 
such a constrained and simple ansatz, we consider the close agreement of the CKM angles and 
quark mass ratios with experiment to be successful. Furthermore, it has not been necessary to 
introduce any large artificial hierarchies in our ansatz (which would require similar hierarchies in 
the corresponding moduli fields), to achieve mass eigenvalues and mixing angles close to those of 
the Standard Model. 




CU 8 0.6 0.4 0.2 

V 

Figure 3: sin# 12 as a function of e and rj. 



13 




0.8 0.6 0.4 0.2 

V 

Figure 4: sin 13 as a function of e and 77. 




0.8 0.6 0.4 0.2 

V 

Figure 5: sin 6*23 as a function of e and rj. 



Mass Ratio 


Blue Region 


Red Region 


Experimental Value 


m c 
m t 


2.5 x 10" 3 


5 x 10~ 3 


5.57- 8.27 x 10" 3 


m 3 
m h 


0.029 


0.041 


0.018 -0.0387 


m u 
m r 


2.4 x 10~ 2 


2.1 x 10~ 2 


1.07-4.5 x 10" 3 



Table 1: Quark mass ratios 
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Figure 6: 5/n as a function of e and rj. 



0.42 



. 4 



\ 



rj 



. 31 



0.36 



0.34 



0.32 



0.3 



.09 



0.1 



0.11 

e 



. 12 



. 13 



Figure 7: The shaded regions highlight the values of e and rj which give an agreement with 
experimental data. The blue/red region without/with taking into account RGE effects. 
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3.2.3 RGE effects 



The CKM matrix is expected to be effected by renormalization group running. For instance, when 
we are close to the so-called quasifixed point, the top Yukawa is large at the GUT or string scale 
and runs to lower fixed point values, and hence RGE effects are important. 

Angle/ Angle GUT 
1.1 r 



0.9 



0.8 



0.7 




Log (M x /Q) 



Figure 8: CKM running from M 



x 



scale. Where red = 



sin 023 



(sin #23 )gut 



blue = 



)16/ 

sm 6 13 
sin9i2 



2 x 10 lb GeV to M w for a top Yukawa of \ t = 5 at the GUT 

sin 9 12 does not run significantly. 



/ sin »i3 



>reen = 



However, the effect on the CKM angles is generally small, as can be seen by inspection of 
the RG equation. Figure || shows how the numerical values of the CKM parameters change with 
renormalization scale when the top quark Yukawa is close to the quasifixed point. The CKM 
angles are suppressed by the running and from the figure we see that the maximum reasonable 
suppression of the 9 2 3 and 9 13 angles is ~ 30%. The phase 8 does not change significantly and is 
expected to be close to |. 



4 Summary and discussion 

We have presented a bottom-up supersymmetric D-brane model with phenomenologically viable 
CP violation, broken by discrete torsion. Furthermore, a simple assumption about the form of 
the kahler metric, motivated by a choice of discrete torsion, produces a CKM matrix described 
by only two free parameters. As a consequence, we predict a single CKM mixing angle and the 
CKM phase to be close to tt/3, both of which lie within current experimental limits. We believe 
this to be the first such prediction from any string model. 

We believe that this class of models may be a first step towards a solution of the SUSY flavour 
and CP problems. Generally, these problems arise because supersymmetry breaking "knows" 
about CP and flavour. The approach that we favour is therefore to use the 'bottom-up' approach 
to build a supersymmetric MSSM, complete with all flavour and CP structure, before discussing 
supersymmetry breaking. 
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An interesting hint may lie in the fact that the Yukawa structure in these models seems to 
favour a hermitian flavour structure. Hermiticity has been suggested as a way to combat large 
EDMs and solve the SUSY CP problem [31[| . In our analysis only the up-quark Yukawas are 
hermitian and so we cannot yet claim to have a solution. However we think that the appearance 
of hermitian flavour structure is intriguing. 

It would be interesting to investigate the dependence of these results on our choice of discrete 
torsion (i.e. our choice of n), and the resulting ansatz for the kahler metrics. Finally an analysis 
of the possibilities for a phenomenologically viable breaking of supersymmetry has been well 
motivated. All of these questions will be the subject of future work. 
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Appendix 



A Calculation of Chan-Paton factors 

Consider a generic chiral superfield whose Chan-Paton factor, as a result of the 70 projection, is 
of size sn^ x suP'. If we act on this state with an element g = (a,b) = uj^uj^ e Z M x Z M , we 
require f\ 

X %xsnU) = r(fc) ( fl ' b h9,3^ (i)xsnU) ll3, (63) 

where r^ k \a, b) is a phase from the action on the string worldsheet fields and for a string endpoint 
on a set of sn^ D-branes, 

i 9 = ©Uv^r a VUJ a Kh. 2 ) 6 ® /(.), (64) 



Im 



where % = e 27 "/ M , Yli m n im = n ^ an d l,m — 0, (M/s) — 1. It follows that we can factor the 
Chan-Paton matrix by writing, 



A inW xsn(i) _ ©*,m^fxl ® Y% (j) . (65) 

lm Im 



where is determined by the action of Z M xZ M and V'*' is left undetermined. Substituting flS5[) 
into (]63|) removes the dependence on giving, 

X (k) tf 1 ti a = r< k Ha,b)fayt a Xto. (66) 

Since the 7^ matrices form a projective representation and hence satisfy 7 9 7h = f3(g,h) r jh'yg, we 
can solve fl66|) with the ansatz X^ k > = 7^7^ ■ This results in, 

/3(u 1 ,u 2 ) Pk ^ 9ka = ^ k \a,b). (67) 



6 This is just a more general version of the projection equations. 
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For 33 sector fields we have, 

r ( £(a,b) = e 2ni{av ^ +bv ^\ (68) 

and for 37 sector fields we have, 

r { $(a,b) =e™( a ^*< +6E ^2). (69) 
Substituting these expressions into fl67j) we find, 

& For 33 fields 



Pk 



(lk 



i£^ fc 4 For 37 fields ' 
For 33 fields 



(70) 



where v^. = jj(a}, af, af). Finally we have, for the 33 fields, 

'(1) _ "n V(2) _ „ 

(3) 



For 37 fields ' 

x^ = 7 -», xf = 7 i l7 |, xf = 7 ~», (7i) 

and for the 37 fields, 

x 37l = 7 jh x 372 = 7 j^ 7 j 2 ^, x 373 = 7 Jr. ( 72 ) 
B Determination of zeros in K33 

As discussed in the main text, the zeros in K33 are inherited from the hermitian matrix, 

Zy = Tr((xg)tX#). (73) 

However, equation (|7lD does not determine the X33 uniquely. Therefore, we choose the n th 
principal root of / y UJk defined by, 

74 = ^(^7)^, (74) 
with S the unitary matrix which diagonalises 7wjt and, 

« a9 h = ^\f. (75) 



Here Aj is the i th eigenvalue of 7aJfc , A^ n = e™ log( ^ A ^ and we restrict arg(Aj) to the interval (— tt, 
thus taking the principal value of log(Aj). 



With these choices, the X33 are uniquely determined and we have, 

^12 - - s \.l^i=i e +e 2^= 9 ( s )+i e ) / 7fi x 

v /V^P™:. 2-rrip/n -27ri(.j-l)/M\ V ' u / 
A lZ^p=0 e Z-/«(p)>j>J(p) e J' 

2iri(l—l)/ns i Jlmln ST^ S „—27ri(l—l)/ns\ 

^13- - s \l^l=i e + e l^l=g(s)+i e ) ( 71 ^ 

v /V^Pmaa: -2mp/n \^ 2ni{j-l)/M\ \ ' ' 1 
A lZ^p=0 e l^u(p)>j>l(p) c !i 

18 



and 




l=g(s)+l c 
Ani(j-1)/M\ 



(78) 



A lZ^ P =o e Z^u( P )>j>i( P ) e 



where, 



Pmax 



r(s— l)n _|_ 11 
L M ^ pi 



u{p) = 

Kp) = 
g(s) = 



(2p + l)f + 1, 



(79) 



fc if s=2k, 



+ 1 ifs=2k+l, 



with [..] denoting rounding to the nearest integer. Furthermore, since the X33 are unitary s x s 
matrices, we have 



Finally, it can be easily seen from (|76|), ([77]) and (f78| ) that Zi 2 = Z i3 = Z23 = for n = 1, 
Z 12 = Z 2 3 = and Z 13 7^ for n = 2 and for n > 3 all elements of Z are non-zero. Hence giving 
rise to the three distinct cases for the distribution of zeros in K 33 as described in the main text. 
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